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WAVE PROPAGATION AT OBLIQUE SHOCKS: HOW DID TYCHO GET ITS STRIPES? 
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ABSTRACT 

We describe a new model for the “stripes” of synchrotron radiation seen in the remnant of Tycho’s 
supernova. In our picture, cosmic rays streaming ahead of the forward shock generate parallel propa¬ 
gating (with respect to the local magnetic field direction) circularly polarized Alfven waves that are 
almost free of dissipation, and due to being circularly polarized exhibit no spatial variation of magnetic 
held strength. Following interaction with the SNR shock with nonzero obliquity, these parallel prop¬ 
agating waves become obliquely propagating, due the the wave refraction (different in principle for 
the different plane wave components), and dissipation sets in. The magnetosonic polarization decays 
faster, due to transit time damping, leaving only the Alfven mode. This surviving mode now exhibits 
a spatial variation of the magnetic held, leading to local maxima and minima in the synchrotron 
emission, i.e. the stripes. We attribute the initial wave generation to the Bell instability, which in 
contrast to the resonant generation of upstream Alfven waves, gives rise to a preferred wavelength, 
and hence the single wave period at which the stripes are seen. Based on estimates for damping rates 
due to turbulent cascade and transit time damping, we estimate the dependence of the visibility of 
the stripes on the shock obliquity, and determine a maximum cosmic ray energy in Tycho’s SNR in 
the range 6 x — lx 10^^ eV. 

Subject headings: acceleration of particles — cosmic rays — magnetic helds — shock waves — ISM: 
supernova remnants 


1. INTRODUCTION 

The advent of the Chandra and XMM-Newton X-ray 
astronomy missions has revitalized the field of cosmic ray 
acceleration. Chandra’s combination of arcsecond reso¬ 
lution imaging and the moderate energy spectral resolu¬ 
tion afforded by CCD X-ray detectors has proven well 
suited to the study of supernova remnants, and specifi¬ 
cally the X-ray synchrotron radiation emitted by cosmic 
ray electrons. This often reveals itself as thin rims of 
continuum emission in the 4-6 keV waveband (a re¬ 
gion essentially free of quasi-thermal line emission from 
shocked plasma) spatially coinciding with the supernova 
remnant forward sho ck. This can be understood (e.g. 

I Vink fc Lamind [200311 as cosmic ray electrons radiating 
in the strong (i.e. amplified) magnetic field at the shock. 
The radial extent is limited by either the radiative loss 
time of the electrons, or the decay of magnetic field post¬ 
shock. The remnant of Tycho’s supernova shows an even 
more intriguing pattern o f synchrotron radia tion; a se¬ 
ries of “stripes” or ripples, (|Eriksen et al Mill), observed 
most clearly in the 4-6 keV waveband where line emis¬ 
sion from the quasi-thermal plasma in the SNR is absent. 
The wavelength corresponds to the gyroradius of cosmic 
ray protons with energies in the range 10^"^ —10^^ eV, but 
the precise mechanism by which these structures form is 
not known. 

The Bell instability (| Belli 12004 120051) , which ampli¬ 
fies magnetic field in the shock precursor region, gener¬ 
ates linearly polarized structures in a near perpendic¬ 
ular shock geometry. Considering the effect of these 
short wavelength fluctuation s on the cosmic ray cur- 
rent, IVladimirov et al.l (|2009D and iBvkov et al.l (|2011al) 
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have shown that long wavelength upstream structures 
can result, with spatial variations of the magnetic field 
strengt h and hence synchrotron emissivity. iBvkov et al.l 
(IMibl) argue that these long wavelength structures 
are responsible for the “stripes”. A number of condi¬ 
tions must be met. Most importantly, the shock re¬ 
gion whe re the stripes appea r must be ‘‘‘nearly perpen¬ 
dicular' (jBvkov et al.ll20lfbl don’t specify how close to 
90° they require), and that in this nearly perpendicu¬ 
lar region, shock acceler ation must still b e efficient. But 
as di scussed elsewhere (iZank et al.l 120061: iLaming et al.l 
1201,?!) . the efficiency of shock acceleration at quasi- 
p erpendicular s h ocks i s open to question. 

iMalkov et al.l (j2012l ) offer an alternative idea, that is 
appears most promising at parallel shocks. They find 
fully nonlinear exact ideal MHD solutions supported by 
the cosmic ray return current, which comprise pulses 
of Alfve n waves that c an pro pagate ahead of the main 
shock. IMalkov et al.l (|2ni2f l argue that these Alfven 
pulses when visible in X-rays should appear as quasi- 
periodic stripes, with spacing simil ar to that of the ob - 
served stripes. In their equations, IMalkov et al.l (I2012D 
neglect thermal and cosmic ray pressure gradients, as¬ 
suming that the ponderomotive force of the turbulence 
is much stronger. In the opposite limit high frequency 
sound waves are g e nerated by the Drury insta bility 
(jPrurv fc FalldllQ^ . iCaprioli fc Spitkovskvl (I2013D sug¬ 
gest that the forward shock may push forwards in up¬ 
stream cavities created at the sa turation of the Bell in - 
stability, also at parallel shocks. iR.akowski et al.l ()2nilD 
discuss a similar idea in connection with shock structures 
(although not stripes) seen in SN 1006. 

In this paper we pursue a different model for these 
synchrotron stripes in Tycho’s SNR. Cosmic rays drift¬ 
ing ahead of a quasi-parallel shock generate upstream 
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Fig. 1. — Schematic diagram of wave propagation at the oblique 
shock, in the shock rest frame. Upstream (left), plasma moves to¬ 
wards the shock front with velocity ui, carrying with it magnetic 
field Bi at angle 9i to ui (or the shock normal). Postshock, the 
magnetic field is B 2 at angle 62 to the shock normal and the flow 
velocity is U 2 at angle a to the shock normal. In the preshock 
medium, waves are assumed to be parallel propagating along Bi. 
Postshock, the waves refract and travel at an angle 92 w th® 
shock normal. They are no longer parallel propagating. In princi¬ 
ple, Alfven and magnetosonic polarizations will refract at slightly 
different angles. The difference between these two angles of re¬ 
fraction tends to zero as the Alfven Mach number 00 , the 

approximation taken in this paper. 

circularly polarized Alfven waves. We assume these 
to be parallel propaKating ( see e.g. iBamert et al.ll200l : 
iGargate fc SDitkovskvll201^ . and with no spatial varia¬ 
tion in magnetic field pressure. Upon passage through 
the shock, both the magnetic field and the wave prop¬ 
agation change direction, and by different amounts, so 
that the formerly parallel propagating waves are now 
obliquely propagating. The circular polarization decom¬ 
poses into its linearly polarized constituents. The mag¬ 
netosonic polarization becomes compressive, and quickly 
decays by transit time damping. The Alfven polarization 
survives longer, and its spatially varying magnetic field 
gives rise the regions of enhanced synchro tron emission, 
observ able as the stripes or ripples seen bv lEriksen et al.l 
(j2011tl . The following sections treat the wave refraction 
at the shock, the wave transmission and reflection coeffi¬ 
cients, and the postshock damping of the magnetosonic 
and Alfven modes. Section 5 puts these topics together 
to explain the origin of the stripes, and to derive physical 
implications from this identification. 

2. WAVE REFRACTION AT SHOCKS 

iMcKenzie fc Westohall (|1970D and I Webb et~aI1 (|1999(1 
give formal accounts of wave properties at oblique 
shocks. We follow and extend slightly for our par- 
ticular application the niore p edagogical treatment of 
lAchterberg fc BlandfordI (jl98(Tll . A schematic diagram 
of the wave propagation at an oblique shock is shown in 
Fig [T] In the preshock medium, we take both Alfven 
and magnetosonic polarizations to be parallel propa¬ 
gating (with respect to the upstream magnetic field) 
with phase velocity Vai, the upstream Alfven speed. 
Downstream, the Alfven polarization has phase speed 
Va 2 cos (02 — (^ 2 w)y where 02 is the angle from the shock 
normal to the postshock magnetic field, and 02 w is the 
angle from the shock normal to the wave propagation 
direction. The phase speed of the magnetosonic wave 
depends on the plasma beta (the ratio of gas pressure to 


magnetic pressure), and is given by fe.g lMelroselll986ll 



+ Cs2 


± 


(^42 + Cs2)^ - 4Ui2C«2 COS^ (02 - 02w) 
Va 2 COS^ (02 - 02 w) - [ COS^ (02 - 02 w) 

^s2 

+ cos"^ (02 - 02 w) ] + ■ . ■ /3 >> 1 

- 4(42 + Cs 2 sin^ (02 - 02 w) /3 < 1- 


( 1 ) 

( 2 ) 


In the following we will only consider the leading terms 
in equations 1 and 2, which will lead to the same angle 
of refraction 02w for both polarizations. In principle, 
however, the different polarizations can refract at slightly 
different angles. At the shock transition, we demand 
that the phases of upstream and downstream waves are 
equal, i.e. ojit — ki • r = uj 2 t — k 2 • r, and that the 
components of wavevectors perpendicular to the shock 
velocity are continuous. For magnetosonic waves at low 
P these constraints give the equations 


^1 Vai - fci cos 0iUi = /c 2 Va 2 - ^ 2^2 cos (02w - a) (3) 
kism0i = fc 2 sin 021V, (4) 

where tana = U 2 ±/u 2 \\ = r’M 2 ±/wi||, the tangent of the 
angle the postshock flow makes to the shock normal with 
r being the shock compression ratio. Dividing (3) by (4) 
to eliminate ki and ^2 leads to 

Vai , „ Va2 , ^ 

---Ml cot 01 =--- U 2 cot 02 W cos a — U 2 sin a, 

sin 01 sm 02 W 

(5) 

which can be rearranged to give an equation for 02 W] 


For Mai —>■ 00 , Ma 2 —>■ 00 , (and hence tana —>■ 0) and 
well away from the precisely perpendicular shock where 
cos 021 V = cos 01 = 0, this reduces to 


cot 021 V = r cot 01 =-— 

tan 02 


(7) 


where we have also used the result tan 02 = r tan 0i from 
equation A4. Writing tan 02iv = tan (O 2 — A) where A 
is the angle between the magnetosonic wave propagation 
direction and the postshock magnetic field direction, we 
find 


. - 1 ) 

tan A = - - cos 0i sin 0,. 

r 


( 8 ) 


The maximum deviation of the wave propagation from 
the magnetic field, and hence the maximum wave com¬ 
pression is achieved for 0i = 45°, where tan A = —15/8 
(assuming r = 4) and A = 62°. Solving for 02 we find 
02 = 76° and 02iv = 14°. 

For magnetosonic waves at high /3 and for Alfven waves 
equation 3 above becomes 


klVAl - fcl cos 01 Ml = A:2 Va2 cos (021V - 02) 
— k 2 U 2 cos (021V — a) 


(9) 
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accounting for the different dispersion properties of the 
Alfven wave. Following the same approach as above, we 
find 


Ma 2 - COS 02 

Ma2 
— tan a + 


cot 92 w = f cot 9i — 


Mai sin 01 


sin 92 
Ma2 ’ 

which can be recast as 
cos 01 


( 10 ) 


1 - 


— tan a + 


Mai 

sin 01 
Mai 


cot 02 W = 'c cot 01 — 


Mai sin 01 


( 11 ) 


In the limit Mai —t oo, ensuring (3 » 1 in the postshock 
medium, this reduces to cot 02 vv = t cot 0i as above. 
Away from this limit, cot 02rv will be different for magne- 
tosonic and Alfven modes. In supernova remnants, this 
is unlikely to be a significant effect, but in the lower Mai 
shocks driven by solar coronal mass ejections, this could 
be an important consideration. The limit Mai oo 
also captures the case when wi —>■ 0, appropriate for the 
Bell nonresonant instability where modes grow at zero 
frequency. 


3. WAVE TRANSMISSION AND REFLECTION AT SHOCKS 

IVainio fc Schlickeised (|1999f) give transmission and re¬ 
flection coefficients for forward and backward propagat¬ 
ing Alfven waves at a parallel shock. Using the jump 
conditions for the tangential electric field, 

\ujiI3i — 0, 

the continuity of the transverse momentum, 

[pUnUt - BnBt/^TT] = 0, 
and the continuity of the mass flux 

[pUn] = 0, 

their result is 

T ^ (MaI+H) (y/fffiljy/F 
R 2 {Mai ± H^) 

Here Ma is the Alfven Mach number of the shock, and 
H is the cross helicity of the waves. iJ = I for for¬ 
ward and H = —1 for backward propagating waves. At 
a parallel shock, the two polarizations, Alfven and mag- 
netosonic, behave the same. At an oblique shock, dif¬ 
ferences emerge. For the Alfven mode, which perturbs 
magnetic field and velocity in the direction perpendic¬ 
ular to the plane in which the wave refraction occurs, 
results for T and R are obtained from equation 15 with 
the replacement Mai Mai / cos 0i. 

The magnetosonic mode is much more involved since 
perturbed magnetic field and velocity vectors lie in the 
plane of refraction. It appears to be tractable only in the 
limit Mai —t oo, when the motion of the shock front itself 
in response to the passing wave can be neglected. We 
evaluate the jump conditions for the tangential electric 
field and the continuity of transverse momentum using 

= — i?2 COS02 -I- ^i?sin02iy (16) 


( 12 ) 

(13) 

(14) 

(15) 


Bt = B2 Su\92 + 5B cos92w (17) 

Un = —W2|| + SuSUi92W (18) 

Ut = U 2 ± + 5u cos 92W ■ (19) 

Substituting into the first jump condition gives 
Ml + 6 u{ sin0i^ (^Bi sin01 + SB( cos0i^ 

— (^—Bi cos 01 -I- 6 B( sin 0i^ Su{ cos 0i (20) 

= (—11211 + 3 u 2 sin02iv) {B2 sin 02 + SBs cos02iv) 

- {-B 2 cos 02 -I- 6 Bs sin 02w) (m 2 ± + Sug cos 02w) • 


Multiplying out and assuming —uiBism9i 
—U 2 \\B 2 sin 02 — U 2 ±B 2 cos 02 yields 


SuiBi — uiSBi cos 01 = SU 2 B 2 cos (02 — 02w) 

— M2||l5il2 COS02W — U 2 ±SB 2 sin02w- (21) 

We put 6ui = 6u{ or Su\ to represent an initially forward 
or backward propagating wave, and then put 6 u 2 = 6 ul+ 
5u\ = ( —SbI + SB 2 ] I\/ 4 :TTp 2 and rearrange to find 


Mai cos 0i ± I = | cos (02 - 02w) 

A" ( ^JrBl 


^A1 „ ,, U2± . 

-cos 02W + Ma\ -sm 02W 

r Ml 

R^ f B 2 f n \ ^A 1 a 

rrh \ <^0® (^2 - d2w) -cos 02W 

i “ ( ^JrBi r 


U2± . „ 

— Mai -sm 92 w 

Ml 


( 22 ) 


The -bve sign in the left hand side corresponds to T/ = 
SB^/SbI and = SB^/SbI, while the -ve sign corre- 
spends to R^ = dB^/SB^ and = 5B\/5B\. A similar 
procedure for the second jump condition gives 


Tf {Mai 

± Myii cos 01-I-cos 201 = < -cos 021V 

R° { yjr 


-MAi\/r— sin 0211 / + ^ cos (02 + 92 w) \ 

Ml Bi j 


Rf 

J'b 

B, 


— p- cos 021 V + Mai sin 02 iv 

Vr Ml 


cos (02 -I- 021V) 


(23) 


Taking Mai —t 00 , so that also M 2 ± —>■ 0 in equations 
22 and 23 we find 


cos 01 


■Vr{Vr±l) 


R^f 2 cos 021 V 

\/sin2 0i+r2 cos2 0i, 


(24) 


where we have substituted f rom cot 02 iv = r cot 0i in the 
final step. This agrees with IVainio fc Schlickeiseil (|I999(1 
and equation 15 in the appropriate limits (0i -A 0 and 
Mai -t 00 ). 
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The forgoing has treated Alfven and magnetosonic 
waves of nonzero frequency in the upstream and down¬ 
stream shock regions. As mentioned above, when being 
driven by the cosmic ray current, these waves grow in the 
upstream region at zero frequency, and hence 5ui =—>■ 0. 
In equations 22 and 23 this leads to the terms ±1 and 
±Mai cos 01 on the left hand sides of equation 22 and 23 
respectively being dropped. We give the corresponding 
expression for the transmission and reflection coefficients 
for Alfven (A) and magnetosonic (M) modes 

T ^ Mait/ cos 01 Agy/r ^ r 
R~ 2 ( Mai/cos0i Agy? ) 2 ^ ’ 

= ]^\Jsin^ 01 -b 7-2 cos2 01 (M), (25) 

which are independent of the cross-helicity H in the limit 
that Mai —>■ oo. 


4. WAVE DAMPING 

Behind the shock, the newly obliquely propagating 
waves are subject to various damping mechanisms. Both 
Alfven and magnetosonic polarizations decay by turbu¬ 
lent cascade, and the magnetosonic wave, now compres¬ 
sive due to its oblique propagation also decays by transit 
time damping by the shocked quasi-thermal ions. 

Many formulations ex ist for the damp ing by turbu¬ 
lent cascade. We follow iBoldvrevI (l2005li who gives an 
expres sion to cover the cases of both the weak turbu- 
lence (iGoldreich fc Sridha iillQQSD a nd strong turbulence 
(|Iroshnikovl 119631 iKraichnanIIl965f) . In these limits the 
approximate damping rates for the large scale turbulence 
are 


7GS = fc2J.^U2 = w tan (02 — 02w) Su2IVa2 (26) 
liK = k 2 ±Su 2 /VA 2 = w tan (02 — 02w) SU 2 IVX 2 , (27) 

where uj = k 2 \\VA 2 from equation 1. The strong turbu¬ 
lence expression for ■jjx is probably the most applicable. 
These simplest expressions refer to balanced turbulence, 
with equal wave amplitudes propagating in each direc¬ 
tion. In our case, and decay rate of the transmitted waves 
will depend on the intensity of counter propagating re¬ 
flected waves, and vice versa. We will assume cosmic rays 
streaming ahead of the shock generate wave travelling in 
one direction only (away from the shock), and then we 
will take du 2 = SuiR, where R is the reflection coefficient 
calculated above in equation 24, in the damping rates. 

The magnetosonic polarization is also subject to tran¬ 
sit time damping by shocked but quasi-thermal ions. 
Thi s damping rat e in the small gyroradius limit is given 
by (IMe]rosell198^ 





A./2sin0i 


'' 


Fig. 2.— Schematic diagram of a plane parallel wave front en¬ 
countering the shock. Upstream the wave is a parallel propagating 
circularly polarized wave, with no spatial variation in the mag¬ 
netic field strength. Behind the shock the magnetosonic polariza¬ 
tion is damped quickly, leaving a spatial variation in magnetic field 
strength due to the surviving Alfven polarization, with maxima ev¬ 
ery half wavelength along the wave propagation direction, Ai/2, as 
indicated. When projected onto the shock front the distance be¬ 
tween magnetic field maxima becomes Ai/2sin0i, where 9i is the 
shock obliquity, the angle of the upstream unperturbed magnetic 
field to the shock normal. 


magnetic field. We represent the ion distribution func- 
tion / as a “kappa” distribution; 


r ^ n _ r(K) r _J^_\ 

(271)^^^^ ^ (k - 3/2) k3/2 [ 2K^VTh2\ 

(29) 

where m is the ion mass, VTh 2 is its thermal velocity, 
r (/c) is the Gamma function of with argument k, and 
K is the index of the distribution. For k —>■ c», / tends 
to a Maxwellian. For finite k, the kappa distribution 
has extended wings compared to a Maxwellian, and may 
be taken to represent suprathermal ions in addition to 
the Maxwellian core. Normalizing the kappa distribution 
over an infinite momentum range, we require k > 3/2 to 
keep the number density n finite, and k > 5/2 to keep 
the energy finite. Smaller values of kappa require a high- 
momentum cutoff to keep particle number and energy 
densities finite. Substituting equation 29 into equation 
28 we get 


277^(7^1/42 1 ^ ’i-r ( k ) 

^11 V TT F (k — 3/2) 

X {2Km‘^VTh2y^^ (30) 



2Km^v^f^2 


k^ 


+ pI 


— K,— l 


p]_dp±. 


ITTD 



fc||U||) 


df 

X hk\i —2Trp±dp±dpi\ , 


(28) 


where uj and are the wave angular frequency and ther¬ 
mal ion gyrofrequency respectively, Va 2 is the postshock 
Alfven speed, and v and p are the ion velocity and mo¬ 
mentum respectively, with subscripts T and || indicat¬ 
ing components perpendicular or parallel to the ambient 


The integral is performed by substituting 
Pj_ = + 2 Km?v^f ^2 tan'0 and evaluat¬ 
ing the resulting integral, sin® t/ cos^'^”® ipdip = 

(1 — cos^ 7 /)^ cos^”"® tjjd (cos^) to find 


7TTD = 


2 TT‘^q^nVl 2 k‘]_VTh2 

[2 

r 1/2 1 

1 1 

fc|| \ 

1 TT 

2 Kv‘^h 2 . 


r(«) 


8M 


'r(K-3/2)K3/2 (k- 1) (k-2)' 


(31) 
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Fig. 3. — Variation with shock obliquity 9bn of the intensity 
contrast A7 between stripe maxima and minima, for the cases of 
Iroshnikov/Kraichnan (IK) or Goldreich-Sridhar (GS) turbulence. 

As K —>■ oo, 7 ttd —t exp (—1/^) /fc||. For 

application below, where /3 >> 1 and w ~ fc||VA2, we 
rewrite as 

ITTD — tan^ (02 - S2w) VTh2lVA2- (32) 

For a fast shock neglecting energy losses to cosmic rays, 
VTh 2 — Vivs/A ^ 2000 km s“^, where the shock velocity 
in Tycho Vs — 5000 km s“^. 


5. HOW DID TYCHO GET ITS STRIPES? 

lEriksen et all (1201111 observe “stripes” in synchrotron 
emission from Tycho’s SNR. The stripes are separated by 
about 10”, which corres ponds to 6 x 10^'^ cm at the pre¬ 
sumed 4.0 kpc distance (iHayato et al.ll201^ . or 5 x 10^^ 
cm at 3.2 kpc (iSlane et al.l 1201411 . We argue that these 
structures arise at oblique shocks, where upstream par¬ 
allel circularly polarized waves which are undamped be¬ 
come obliquely propagating damped waves postshock. 
The key point is that the magnetosonic polarization 
damps faster than the Alfven polarization, and so a spa¬ 
tial variation in the magnetic field strength will emerge 
associated with the surviving Alfven mode. Along the 
direction of propagation the absolute magnitude of the 
magnetic field strength has a maximum every half wave¬ 
length (Ai/2). When projected onto the shock front, 
the maxima are separated by a distance of Ai/2sin0i, 
as shown schematically in Fig[^ Although less damped 
than the magnetosonic polarization, the Alfven modes 
still damp within a wavelength or so postshock. These 
magnetic field variations are illuminated by synchrotron 
radiation from cosmic ray electrons, the radiative cooling 
time for which typically lies between the magnetosonic 
and Alfven mode damping times. This hierarchy en¬ 
sures that the damped magnetosonic polarization does 
not contribute significantly to the emission, while the 
Alfven mode only contributes close to the shock, while 
the electrons are still radiating in X-rays. 

In the case of Iroshnikov-Kraichnan turbulence, the 
difference in synchrotron intensity between maxima and 
minima in the stripes may be written 


A/oc 


exp - exp {-'JIKt - 

X cos'’+^ 01 


(33) 


where s is the index of the synchrotron photon spec¬ 


trum and s -f 1 is the power law dependence of the mag¬ 
netic field strength on the synchrotron emission. The 
final factor of cos0i gives the approximate shock obliq¬ 
uity dependence of the preshock magnetic field amplifica¬ 
tion. Goldreich-Sridhar turbulence simply requires jcs 
in place of 'Jik- This intensity contrast is maximized at 
t — In (1 -I- 'yrTo/liK) IjttdI (s + 1) with value 


AI oc 


exp {—jik/jttd) In (1 -I- 'yTTo/liK) 


— exp {—jik/jttd — 1) In (1 -|- jttd/jik) 
xcos"+i0i. (34) 


Figure 3 shows the dependence of AI on the shock obliq¬ 
uity 01, for both Iroshnikov-Kraichnan (IK; strong) tur¬ 
bulence and Goldreich-Sridhar (GS; weak) turbulence, 
calculated from equation 34, using equations 26, 27, and 
32. We take Su 2 = SuiR, and Sui ~ 180 km s“^, cor¬ 
responding to i? ~ 2 from equations 25, and postshock 
magnetic fiel d and density of 1 80 /rG and 1.2 amu cm“^ 
respectively (jSlane et al.ll20l4) . In both cases the emis¬ 
sion is relatively strongly peaked over a restricted range 
of angles, possibly suggesting a reason why the stripes 
are not ubiquitous over Tycho, but only seen in certain 
“special” regions. 

The separation at which the stripes preferably appear 
can be used to estimate the maximum cosmic ray ion 
energy in Tycho. The Bell instability grows structures 
at a para llel wavevector i n the ambient mag netic field 
given by (lBellll200l . 120051 iLaming et al.ll201^ 


fell = 1.5 X 10“® 


rjm 


Ui 


-1 


,5000 km s 

7maa;- 371/4 _ ^ 

7m.aa:7l {\^^max 1) Ai 



(35) 


Here rj is the fraction of shock energy going into cos¬ 
mic rays, nt is the ion number density in the preshock 
medium. The highest energy cosmic rays have Lorentz 
factor jmax, while the highest energy magnetized cosmic 
rays have Lorentz factor 71. The cosmic ray current with 
particles with 71 < 7 < 'jmax drives the B ell instability. 
We ha ve also set 5B/B = 1 in equation 3 of iLaming et al.l 
(|2014f) for feyrg, appropriate for the far upstream region 
of the cosmic ray precursor, and divided through by the 
cosmic ray gyroradius. 

We identify Ai / 2 sin0i = 5 x 10^^ cm at a distance to 
Tycho of 3.2 kpc (jSlane et al.ll20l4 model A), and tak¬ 
ing 01 from Fig 3, we estimate ^max- Letting 71 "(max, 
and taking 77 = 0.26, rii = 0.3cm“^, we find "(max = 
6 X 10® — 1 X 10® corresponding to 0i = 10° — 15°. Al¬ 
though mo del dependent, this broadly agrees with the 
modeling of iSlane et al.l (I2014I1 . who find the cosmic ray 
energy spectrum breaking from something close to a 
power law to a steeper fall-off just a bove 7 ^ 10® 
(their Fig. 4), and the original estimate of lEriksen et al.l 
(j2011ll . who simply identified the spacing of the stripes 
with a cosmic ray gyroradius. Projected onto the shock 
velocity vector, the wave phase varies with distance 
Ai/2 X COS01 ~ 1 X 10^^ — 1.5 X 10^^ cm. Thus with a 
shock velocity of 5000 km s“^, changes in the projected 
position of the stripes should be visible in 2 x 10® —3 x 10® 
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seconds, i.e. 6-10 years. lEriksen et al.l (|2011li report an 
intensity contrast of a factor of 25 between stripe max¬ 
ima and minima. Relating this to with s ^ 2.11 

we find a variation of B between maxima and minima 
of about a factor 3. This is significantly lower than the 
likely upstream magnetic field contrast of between a fac¬ 
tor of 9 and 15, but probably reasonable when we con¬ 
sider that both polarizations are amplified upstream by 
this amount, and then the magnetosonic polarization is 
preferentially damped by transit time damping. In the 
case that the electron cooling time becomes longer than 
the Alfven mode damping time, the variation of magnetic 
field along the shock front will be less clear, but stripes 
separated by A 2/2 oriented normal to k 2 extending fur¬ 
ther inside the shock would result, analogous to ripples 
associated with sound waves. 

6. CONCLUSIONS 

We suggest a model for the sync hrotron “stripes” or 
ripples observed in Tycho’s SNR (|Eriksen et al.l 1201 Ifl 
based only on simple ideas about Alfven propagation and 


dissipation at the forward shock, calculated as a function 
of its obliquity. We find that the contrast is highest, and 
therefore that the stripes would be most visible for a nar¬ 
row range of shock obliquities close to 15° in the case of 
strong Iroshnikov-Kraichnan turbulence. In our model, 
the emission comes from behind the shock front, once 
the magnetosonic component of the originally circularly 
polarized wave has been transit time damped. This is in 
contrast to previous ideas which have placed the struc¬ 
tures giving rise to the synchrotron stripes in the shock 
precursor. In the model advanced here, some motion of 
the stripes over a period of one to a few years should be 
visible, leading to possibilities for an observational test. 
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APPENDIX 

MHD SHOCK RELATIONS 


Here we reproduce a few relationships relating to MHD shock, following mainly iMelrosel (jI986[l . The postshock 
perpendicular flow speed is 

(r-I)sin 6 »iCOs 6 »i 


U2± = ui- 


— r cos^ 9i 


—^ 0 as ^ 00 . 


The transverse magnetic field is 



^ r (Ml, — cos^ 6i) 

B 2 sin O 2 = , ,0 -T-— Bi sin 0i —^ rBi sin 0i as Mai —>• 00 . 

M^i — r cos^ 0i 

(A2) 

Writing 

i? 2 sin 02 r (M^, — cos^ 0i) cos 0, sin 02 

Hi sin 01 M^,—r 008^01 cos 02 sin 0i 

(A3) 

we hnd 

/, (r — I) sin 01 cos 01 ^ 1 

tan 02 = r tan 0i H-p,--— = r tan 0i -I— tan a 

^Ai ~ cos^ 01 r 

(A4) 

to be compared with 

cot 02 VV = r cot 01 — tan a 

(A5) 


from equation 6 , for the wave refraction. 
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